A primer on problems and prospects of dark energy 
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This review on dark energy is intended for a wider audience, beginners as well as experts. It 
contains important notes on various aspects of dark energy and its alternatives. The section on 
Newtonian cosmology followed by heuristic arguments to capture the pressure effects allows us 
to discuss the basic features of physics of cosmic acceleration without actually resorting to the 
framework of general theory of relativity. The brief discussion on observational aspects of dark 
energy is followed by a detailed exposition of underlying features of scalar field dynamic relevant 
to cosmology. The review includes pedagogical presentation of generic features of models of dark 
energy and its possible alternatives. 

PACS numbers: 98.80 Cq 



I. INTRODUCTION 

Twentieth century has witnessed remarkable develop- 
ments in the field of cosmology. The observation of red- 
shift of light emitted by distant objects and the discovery 
of microwave background in 1965 have revolutionized our 
thinking about universe. The hot big bang model then 
received the status of standard model of universe. How- 
ever, in spite of the theoretical and observational suc- 
cesses, cosmology remained confined to a rather narrow 
class of scientist; others considered it as the part of a 
respectable philosophy of science. Cosmology witnessed 
the first revolution in 1980 with the invent of cosmo- 
logical inflation making it acceptable to the larger com- 
munity of physicists. Since then it goes hand in hand 
with high energy physics. The scenario envisages that 
universe has gone through a phase of fast accelerated ex- 
pansion at early epochs. Inflation is a beautiful paradigm 
which can resolve some of the in built inconsistencies of 
the hot big bang model and provides a mechanism for 
generation of primordial fluctuations needed to seed the 
structure we see in the universe today. In the past two 
decades, observations have repeatedly confirmed the pre- 
dictions of inflation. However, its implementation is ad 
hoc and requires support from a fundamental theory of 
high energy physics. As inflation takes place around the 
Planck epoch, the needle of hope points towards string 
theory — a consistent theory of quantum gravity. 

The second revolution cosmology witnessed in 1998, is 
related to late time cosmic acceleration P, The ob- 
servations of high redshift supernovae reveals that uni- 
verse is accelerating at present. The phenomenon is indi- 
rectly supported by data of complimentary nature such 
as CMB, large scale structure, baryon acoustic oscillation 
and weak lensing. It is really interesting that the thermal 
history of our universe is sandwiched between two phases 
of accelerated expansion. In the Newtonian language, 
cosmic repulsion can be realized by supplementing the 
Newtonian force by a repulsive term on phenomenological 
grounds. The rigorous justification of the phenomenon 
can only be provided in the frame work of general the- 
ory of relativity (see Ref.JH for early attempts in this 



direction). Late time acceleration can be fueled either 
by an exotic fluid with large negative pressure dubbed 
dark energy^, 0, 0, @, H, H, G3 or by modifying the grav- 
ity itself [Tl| ■ The simplest candidate of dark energy is 
provided by cosmological constant A[l2|, EH, EH, though, 
there are difficult theoretical issues associate with it. Its 
small numerical value leads to fine tuning problem and 
we do not understand why it becomes important today 
a la coincidence problem. 

Scalar fields provide an interesting alternative to cos- 
mological constant [3, E3]- To this effect, cosmological 
dynamics of a variety of scalar fields has been investi- 
gated in the literature(see reviewQ for details). They 
can mimic cosmological constant like behavior at late 
times and can provide a viable cosmological dynamics 
at early epochs. Scalar field models with generic features 
are capable of alleviating the fine tuning and coincidence 
problems. As for the observation, at present, it is abso- 
lutely consistent with A but at the same time, a large 
number of scalar field models are also permitted. Future 
data should allow to narrow down the class of permissible 
models of dark energy. 

As an alternative to dark energy, the large scale mod- 
ifications of gravity could account for the current accel- 
eration of universe. We know that gravity is modified at 
short distance and there is no guarantee that it would not 
suffer any correction at large scales where it is never ver- 
ified directly. Large scale modifications might arise from 
extra dimensional effects or can be inspired by fundamen- 
tal theories. They can also be motivated by phenomeno- 
logical considerations such as f(R) theories of gravity. 
However, any large scale modification of gravity should 
reconcile with local physics constraints and should have 
potential of being distinguished from cosmological con- 
stant. To the best of our knowledge, all the schemes of 
large scale modification, at present, arc plagued with one 
or the other problems. 

The review is organized as follows: After introduction 
and a brief background, we present cosmology in Newto- 
nian framework in section III and mention efforts to put 
it on the rigorous foundations in the domain of its valid- 
ity. In section IV, we put forward heuristic arguments 
to incorporate A, in particular and pressure corrections, 
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in general, in the evolution equations and describe the 
broad features of cosmological dynamics in presence of 
cosmological constant. In section V, we provide a short 
introduction to relativistic cosmology and discus issues 
associated with cosmological constant. After a brief sub- 
section on observational aspects of cosmic acceleration, 
we proceed to highlight the generic features of scalar field 
dynamics relevant to cosmology and mention the current 
observational status of dynamics of dark energy. In the 
last section before summary, we present a discussion on 
the current problems of alternatives to dark energy. 

Last but not the least, a suggestion for the follow up of 
this review is in order. At present, there exist, a number 
of excellent reviews on dark energy [i, IH 1(1 M. l9l. Hoj and 
cosmological constant [HI, [ill [l4| which focus on different 
aspects of the subject. Four recent and very interesting 
reviews 0, [H, H [H| which try to address the theoretical 
and observational aspects of late time cosmic acceleration 
are highlyrecommended. Humility does not allow to say 
that Ref. Q is the most comprehensive theoretical review 
on dark energy with pedagogical exposition. 



II. THE SMOOTH EXPANDING UNIVERSE 

Universe is clumpy at small scales and consists of very 
rich structure of galaxies, local group of galaxies, clusters 
of galaxies, super clusters and voids. These structures 
typically range from kiloparsecs to 100 megaparsecs. The 
study of large scale structures in the universe shows no 
evidence of new structures at scales larger than 100 mega- 
parsecs. Universe appears smooth at such scales which 
leads to the conclusion that universe is homogeneous and 
isotropic at large scales which serves as one of the fun- 
damental ass ump tions in cosmology known as cosmolog- 
ical principal\l&^. Homogeneity tells us that universe 
looks the same observed from any point whereas isotropy 
means that It looks same in any direction. In gen- 
eral these are two independent requirements. However, 
isotropy at each point is stronger assumption which im- 
plies homogeneity also. Cosmological principal presents 
an idealized picture of universe which allows us to un- 
derstand the background evolution. The departure from 
smoothness can be taken into account through perturba- 
tions around the smooth background. Observations con- 
firm the presence of tiny fluctuations from smoothness 
in the early universe. According to modern cosmology, 
these small perturbations via gravitational instability are 
believed to ha ve g rown into the structures we see today 
in the universe [II [H HH, [S3, [H, [H [H [H, H3 . 

One of the most remarkable discoveries in cosmology 
includes the expansion of universe and its beginning from 
the big bang. The analysis of radiation spectrum emitted 
from distant galaxies shows that wavelengths of spectral 
lines are larger than the actually emitted ones, the phe- 
nomenon is known as redshift of light. The redshift is 
quantified by symbol z defined as, z = (A 6 — A em )/A em . 
According to the Doppler effect, the wavelength of light 



emitted by a source receding from the observer appears 
shifted towards red end of spectrum and the redshift is re- 
lated to the velocity of recession v as, z ~ v/c for v << c. 
In the beginning of the last century, astronomers could 
measure the distances to a number of distant galaxies. 
Hubble carried out investigations of recession velocities 
and plotted them against the distances to galaxies. He 
concluded in 1929 that there is a linear relation between 
recession velocity of the galaxies and the distance to them 
— the so called Hubble Law. 

The observational conclusion that universe expands is 
based upon the redshift of radiation emitted by distant 
galaxies. Can we have another explanation for the red- 
shift? It might look surprising that photons from larger 
distances emitted from galaxies reach us redshifted due 
to the recession of galaxies and nothing else happens to 
them. They travel through intergalactic medium and 
could be absorbed by matter present there and then emit- 
ted loosing part of their energy in this process thereby 
leading to their redshift without resorting to expansion of 
universe. This apprehension can be refuted by a simple 
argument. As for the absorption, the underlying process 
is related to the scattering of photons by the particles 
of intergalactic medium. If true, the source should have 
appeared blurred which is never observed. Other efforts 
assuming the exotic interactions of photons could not 
account for the observed redshift. Thus the only viable 
explanation of the phenomenon is provided by the ex- 
pansion of universe [22j. 

If we imagine moving backward in time, universe was 
smaller in size, the temperature was higher and there was 
an epoch when the universe was vanishingly small with 
infinitely large energy density and temperature-the be- 
ginning of universe dubbed as the big bang. The matter 
was thrown away with tremendous velocity, since then 
the universe is expanding and cooling. At early times it 
was extremely hot and consisted of a hot plasma of el- 
ementary particles, there were no atoms and no nuclei. 
Roughly speaking, at temperatures higher than the bind- 
ing energy of hydrogen atom, the photons were freely 
scattering on electrons and atoms could not form. As 
the universe cools below the temperature characterized 
by the binding energy of hydrogen atom, electrons com- 
bine with protons to form hydrogen atom leading to the 
decoupling of radiation from matter. This was an im- 
portant epoch in the history of universe known as re- 
combination. The decoupled radiation since than is just 
expanding with the expanding universe and cooling. The 
discovery of microwave background, the relic of the big 
bang, in 1965 confirms the hypothesis of hot big bang. 



III. THE HOMOGENEOUS AND ISOTROPIC 
NEWTONIAN COSMOLOGY 

Newtonian theory of gravitation allows us to under- 
stand the expansion of a homogeneous isotropic universe 
in a simple way. Newtonian description is valid pro- 
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vided the matter filling the universe is non-rclativistic 
and scales associated with the problem are much smaller 
than the Hubble radius. For instance, at early epochs, 
the universe was hot dominated by radiation. Hence 
early universe strictly speaking, should be treated by rel- 
ativistic theory. The general theory effects are also cru- 
cial at super Hubble scales. Despite its limitations, New- 
tonian cosmology provides a simple and elegant way of 
understanding the expansion of universe 



A. Hubble law as a consequence of homogeneity 
and isotropy 

Using the Newtonian notions of physics, let us show 
that the Hubble law is a natural consequence of homo- 
geneity and isotropy. Let us choose a coordinate system 
with origin O such that matter is at rest there and let 
us observe the motion of matter around us from this co- 
ordinate system. The velocity field, i.e., the velocity of 
matter at each point p around us at an arbitrary time, 
depends upon the radius vector r and time t. We should 
now look for the most general velocity field in a homo- 
geneous and isotropic universe. Let us assume another 
observer located at point O' with radius vector r > and 
moving with velocity v(r D ') with respect to the observer 
O. If we denote the velocity of point p relative to O and 
O' at time t by v(r p ) and v'(i-p), we have, 

r 'p = T P - r o> (1) 

v'(r;)=v(r p )-v(r ,) (2) 

where r p and r' p denote the radius vectors of point p 
with respect to O and O' respectively. The cosmological 
principal tells that the velocity field should have the same 
functional form at any point, 

v(r;)=v(r p )-v(r o (3) 

which clearly implies that the velocity field is a linear 
function of its argument r, 

v{r,i)=T(t)r (4) 

where T is a 3 x 3 matrix. The matrix can always be 
diagonalized by choosing a suitable coordinate system. 
Isotropy then reduces it to kronecker symbol (Tj j = 
H(t)Sij) leading to 

v(r,t)=JT(f)r (5) 

where H is known as the Hubble parameter. In gen- 
eral, a velocity field can always be decomposed into ro- 
tational part, inhomogeneous part and isotropic part at 
each point. It is not then surprising that the homoge- 
neous and isotropic velocity field has the form {5]) known 
as Hubble law. 



It can easily be verified that Hubble law holds at any 
point. If we move from O to O' , we can write 

v'(r' p ) = Hr p - Hr , = H(t)r' p (6) 

The Hubble's law gives the most general form of veloc- 
ity field permissible by the homogeneity and isotropy of 
space. 

Hubble law tells us how distance between any two 
points in space changes with time provided we know the 
expansion rate given by H(t), 

r(t) = xe-fo H W d \ x = r(t = 0) (7) 

The law of expansion depends upon how the Hubble pa- 
rameter H varies with time. Eq. ([7]) shows how distances 
in a homogeneous and isotropic universe scale with the 
scale factor a(i), 

a(t) = efi or H(t) = - (8) 

a 

r(t) = a(t)x (9) 

The complete information of dynamics of a homogeneous 
and isotropic universe is contained in the scale factor; we, 
thus, need evolution equation to determine a(t). In case 
H is independent of time we have exponentially expand- 
ing universe dubbed de-Sitter space. In what follows 
we shall confirm that constant Hubble rate is allowed 
in relativistic cosmology provided the energy density of 
matter in the universe is constant. It is believed that 
universe has passed through an exponentially expanding 
phase known as inflation at early times. 

According to Hubble law, in a homogenous and 
isotropic universe, all the material particles move away 
radially from the observer located at any point in the uni- 
verse. This motion is refereed to as Hubble flow. Indeed, 
any freely moving particle in such a background would ul- 
timately follow the Hubble flow. Motion over and above 
the Hubble flow is called peculiar motion which can only 
arise in a perturbed universe. It often proves convenient 
to change a coordinate system dubbed comoving which 
expands with expanding universe. Matter which follows 
the Hubble flow will be at rest in the comoving coordi- 
nate system, i.e. matter filling a homogeneous isotropic 
universe is at rest with respect to the comoving observer. 
Both the frames are physically equivalent. Let us clarify 
that universe does not appear homogeneous and isotropic 
to any observer; for instance if an observer is moving with 
a large velocity say towards a particular galaxy, universe 
looks different to him/her. A physical coordinate system 
is the system in which matter is at rest at the origin and 
moves away radially at other points. The radius vector 
r of any point in this system called physical, changes 
with time whereas its counterpart x in the comoving sys- 
tem is constant. This means that physical distance be- 
tween any two points in the expanding universe is given 
by the comoving distance multiplied by a factor which de- 
pends upon time which is precisely expressed by Eq.([7]) 
or equivalcntly by Eq.®. 
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continuity equation to have the usual form, 



FIG. 1: Particle of mass m on the surface of a sphere of radius 
r(i) in an expanding universe with uniform matter density. 



B. Evolution equations 

We now turn to the evolution equation for the scale 
factor. Thank to isotropy, we can employ spherical sym- 
metry to derive the evolution equation. At a given time 
t called the cosmic time, let us consider a sphere cen- 
tered at O with radius r(t). Let pb(t) be the density 
of matter in the homogeneous isotropic space referred to 
as background space hereafter. We assume that the net 
gravitational force on a particle of mass m situated on the 
surface of the sphere due to matter out side the sphere 
is zero which means that matter inside the sphere alone 
can influence the motion of the particle. The total energy 
of the particle on the surface of the sphere (see FigQ]) at 
any time is constant given by the expression [5^], 



1 2 47T 2 

E To t = -mr - —mGp b r 



(10) 



This equation can be cast in the following convenient 
form, 



H 2 
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(ii) 



which readily translates into an evolution equation for 
a(t) (see, Eq.([9])) known as Friedmann equation, 
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(12) 



where K can be zero, negative or positive depending how 
kinetic energy compares with the potential energy. 

In order to solve the evolution equation for a(t), we 
need to know how matter density Pb(t) changes with time, 
i.e., we need the conservation equation in the expanding 
universe. For non-rclativistic fluid, the continuity equa- 
tion that gives us the evolution of matter density of the 
fluid is, 



dp b it) 
dt 



(V.pftv) = 



(13) 



Remembering that the matter density of the background 
fluid is independent of the coordinates and the fluid ve- 
locity is given by the Hubble law ([H]), we transform the 



d Pb (t) 
dt 



3Hp b = 



which formally integrates to, 



P b (t)= Pb {-) 



(14) 



(15) 



where the subscript '0' denotes the quantities at the 
present epoch. The evolution of matter density of non- 
relativistic fluid has a simple meaning that mass of fluid 
in a co- moving volume remains constant. 

Though the Eq. (fT2|) formally resembles the evolution 
equation of relativistic cosmology, its derivation pre- 
sented above is defective. The expression for the potential 
energy is written with an assumption that gravitational 
potential can be chosen zero at infinity which is not true 
in an infinite universe. Since the mass density pb is con- 
stant in space, the total mass of universe diverges as r 3 . 
As a result, the potential — AirGpbr 2 /3 can not be nor- 
malized to zero at r = oo. One could try to circumvent 
the problem by assuming that pb vanishes for a given 
large value of r but it would conflict with the underlying 
assumption of homogeneity. Therefore, conservation of 
energy is difficult to understand in an infinite universe 
with uniform matter density. 

We can also derive the evolution equations using the 
Newtonian force lawpjl . The force on the unite mass sit- 
uated on the surface of homogeneous sphere with radius 
r is given by 



4ttG 
F = pr 



The Euler's equation 
9v 



, v.V)v = ^ + F 

dt p b 



(16) 



(17) 



in a homogeneous isotropic background simplifies to 

¥={H + H 2 )t (18) 

where F is the force per unit mass on the fluid element 
given by Eq. p^|) . We have used the fact that pressure 
gradients are absent in a homogeneous isotropic back- 
ground and the velocity field is given by the Hubble 
law. It should also be noted that the pressure Pb = 
for the non-relativistic background fluid under consid- 
eration. Using expressions p^|) & (fT5)) . we obtain the 
equation for acceleration, 



4ttG , . 

z~na = — r^w 



a dt? 



(19) 



which could also be obtained directly from Eq. (|T6]) . 
Equation (|19[) can easily be integrated to give the Fried- 
mann equation. Indeed, by multiplying the above equa- 
tion by a and using the evolution of mass density allows 
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us to write 



H 1 



K = ai 



InG . . K 
, ( 8irGp 



-m 



(20) 
(21) 



The above derivation is also problematic as it assumes 
that mass out side the sphere, used while writing Eq. (fT6|) . 
can be neglected which is not true for infinite universe 
with constant mass density. 

The problem can be circumvented by using the geo- 
metric reformulation of Newtonian gravity in the lan- 
guage of Cartan. According to Cartan's formulation, or- 
bits of particles are assumed to be the geodesies of an 
affine space and gravity is then described by the curva- 
ture of the affine connectionisee Ref.[30j and references 
therein) . According to Ref. [3(| , no pathology in cosmol- 
ogy associated with Newton's force law then occurs and 
the evolution equations of Newtonian cosmology, 
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8irG . , K 
—z-Pb{t) ~ 



K = al 



1 d 2 a 
a~aW '' 
d Pb (t) 



lirGp 
3 

4irG 



-m 



at 



p b (t) 



3H Ph (t) = 



(22) 
(23) 

(24) 
(25) 



can be put on rigorous foundations. Eqs. (f2"2")) . (j2"4")l . & 
(j25[) are identical to evolution equations of Friedmann 
cosmology for non-relativistic fluid filling the universe. 
Whether or not one adopts the formulation presented 
in Ref. [3(| , Newtonian cosmology is nevertheless elegant 
and simple. 

Let us point out an important feature of Newtonian 
cosmology. We note that the expression of Kjo? remains 
unchanged under the scale transformation a(t) — ► Ca(t), 
C being constant. As a result, the evolution equations 
(|2"4")) & (|2"2"|) also respect the scale invariance. This in- 
variancc is a characteristic of specially flat Friedmann 
cosmology. The Newtonian cosmology can mimic all 
three topologies of relativistic cosmology corresponding 
to K = 0, ±1 in spite of the fact that the underlying ge- 
ometry in Newtonian cosmology is Euclidean. Let us note 
that the scale factor in Newtonian cosmology can always 
be normalized to a convenient value at the present epoch. 
This is related to a simple fact that the Friedmann equa- 
tion (|2"2"1) docs not change if we re-scale the scale factor 
which leaves the normalization of a arbitrary. The often 
used normalization fixes the scale factor a(t) = 1 at the 
present epoch, i.e, ao = 1. In case of relativistic cosmol- 
ogy, the latter can only be done in case of K = whereas 
in case of K = ±1, the numerical value of the scale factor 
ao depends upon the matter content of universe. 

The second important feature of Newtonian cosmol- 
ogy is that it leads to an evolving universe. Indeed, we 
could ask for a static solution given by a(t) and a{t) = 



which is permitted by the Friedmann equation (|22|) but 
not allowed by the equation for acceleration (j2"4")l . It is 
really remarkable that Newtonian cosmology gives rise to 
an evolving universe. It is an irony that the discovery of 
expansion of universe had to wait the general theory of 
relativity. This is related to the commonly held percep- 
tions of static universe which was prevalent before Fried- 
mann discovered the non-static cosmological solution of 
Einstein equations. So much so that Einstein himself did 
not believe in the Friedmann solution in the beginning 
and tried to reconcile his theory with static universe by 
introducing cosmological constant which he later with- 
drew. 



C. The past, the future and how old are we? 

The general features of solutions of evolution equations 
can be understood without actually solving them. What 
can we say about the past and the fate of universe? The 
equation for acceleration tells us that a < for standard 
form of matter. This means that a(t) as a function of 
time is concave downward. We need input regarding a 
at present to make important conclusion about the past. 
Observation tells us that a(t) > at present. Thus a(t) 
monotonously decreases as t runs backward. It is there- 
fore clear that there was an epoch in the history of uni- 
verse when a(i) vanishes identically. Without the loss of 
generality we can take t = corresponding to a(t) = 0. 

As for the fate of universe, the problem is similar to 
that of escape velocity, namely, if K > 0, the kinetic en- 
ergy is less than the potential energy. In this case a(t) 
would increase to a maximum value where a(t) = 0, it 
would start decreasing thereafter till it vanishes and uni- 
verse ends itself in big crunch. In case, K < 0, scale fac- 
tor would go on increasing for ever; K = represents the 
critical case. Three different possibilities, K = 0, K > 
or K < correspond to critical, closed and open universe 
respectively. We should emphasize that the fate of uni- 
verse also crucially depends upon the nature of matter 
filling the universe. In some case, the universe may end 
itself in a singular state or the cosmic doomsday. 

Which of the three possibilities is realized in nature? 
To answer this question, let us rewrite Eq. (|12[) in a con- 
venient form, 



n b (t) - 1 



K 



(aH)' 



:, n t (t) 



Pb(t) 

Pc(t) 



(26) 



where the critical density is defined as, p c (t) = 
3H 2 (t)/8irG. Specializing the expression (|26|) to the 
present epoch, we find that, 



(0) 

b 

(0) 
b 

(0) 



> i (/4 0) > /4 0) 



HpT-p^) 



(0)1 



<i(p! 0) <p^) 



K > — » closed universe, 
K = — > critical universe, 
K < — > open universe. 



where the super script '0' designates the correspond- 
ing physical quantities at the present epoch. Since we 
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know the observed value of pf\ one of the three types 
of universe we live in, depends upon how matter density 
in universe compares with p^ ■ Observations on Cosmic 
Microwave background (CMB) indicate that universe is 
critical to a good accuracy or K ~ which is consistent 
with inflationary paradigm. 

Let us come to the solution of Newtonian cosmology 
in case of K = 0. Substituting Pb(t) from Eq. (fT5l) in 



Ea.(f2"2"]). we find that, a 2 
giving rise to 



which easily integrates 



a(t) = ( - 



2/3 



Pb(t) = P b 

w 3t 



(0) {to 



(27) 



(28) 
(29) 



(30) 



The above solution is known as Einstein-de-Sitter solu- 
tion. We can estimate the age of universe using Eq. ([29|) . 

_ 2 1 

Interestingly, if gravity were absent, universe would ex- 
pand with constant rate given by Hq. Using Hubble law 
we would then find, 

1 

Ho 



in 



(31) 



which is the maximum limit for the age of universe in the 
hot big bang model (2H ~ 1 /3 < t < H- 1 ). The presence 
of standard matter always leads to deceleration thereby 
leading to smaller time taken to reach the present Hubble 
rate of expansion. The presence of cosmological constant 
or any other exotic form of matter can crucially alter this 
conclusion. 



D. Cosmological constant a la Hooke's law 

We have seen that Newtonian cosmology gives rise to 
evolving universe but for the historical reasons, cosmol- 
ogy had to wait the general theory of relativity to dis- 
cover it. The fact that Newtonian cosmology leads to 
non-stationary solution was known before general theory 
was discovered but it could receive attention as it con- 
flicted with perception of static universe. Attempts were 
then made to modify Newtonian gravity to reconcile it 
with the static universe. Clearly, the modification should 
be such that it becomes effective at large scales leaving 
local physics unchanged. Looking at the Newton's force 
law ([16]) . it is not difficult to guess that static solution 
is possible provided that we add a repulsive part propor- 
tional to the radius vector r in Eq. fll)]) , Newton's law of 
gravitation should therefore be supplemented by linear 
force law[H[23,[3ll,[13 



4ttG 1 A 

■— P h r+-Ar 



(32) 



where A is known as cosmological constant which is pos- 
itive in the present context. It is interesting to note 
that there are only two central forces namely, the inverse 
square force and the linear force which give rise to stable 
circular orbits. 

Our discussion of cosmological constant is heuristic 
and the cheap motivation here is to incorporate the re- 
pulsive effect in the evolution equations. We rewrite the 
modified force law ((32)) as an equation of acceleration 
using the comoving coordinates, 



1 d 2 a 
adt^ 



AirG 



Pb(t) 



A 

3 



(33) 



which shows that a positive A term contributes to accel- 
eration as it should. The integrated form of Eq. ([33|) is 
given by, 



H 



8ttG . N K 



A 
~3 



(34) 



where the integration constant K can be formally written 
again through physical quantities defined at the present 
epoch. The modified force law (|32p was proposed much 
before Einstein's general theory of relativity by Neumann 
and Seeliger in 1895-96(11 HI- 

Let us note that adding cosmological constant to New- 
tonian force is equivalent to adding a constant matter 
density p\ = A/8itG to the background matter density 
Pb which does not to go well with the continuity equation 
(Ti"4|) . Since the acceleration equation also gets modified 
in presence of A, we should check whether the modified 
evolution equations allow this possibility. If we differenti- 
ate Eq. (|3"4"|) with respect to time and respect the modified 
acceleration equation, we find that constant matter den- 
sity is permissible in the expanding universe. As for the 
continuity Eq. (jT4j) , it is valid for a perfect non-relativistic 
fluid. The cosmological constant does not belong to this 
category, the pressure corresponding to constant energy 
density is not zero. The continuity equation should take 
the note of pressure and get appropriately modified. As 
pointed out earlier our present discussion of cosmological 
constant here is qualitative. Rigorously speaking, we are 
trying to get the right thing in the wrong place! We shall 
come back to this point after we incorporate the pressure 
corrections in the evolution equations. 

Evolution equations (|34|) and (|33|) admit a static solu- 
tion (a = const = an) in case of K > 0. Static Einstein 
universe ( d = and a = 0) is possible provided that A 
has definite numerical value 



A = K 



4nGp 



(0) 



(35) 



We shall observe after a short while that the static Ein- 
stein universe is unstable under small fluctuations. 

The qualitative features of solutions of evolution equa- 
tions can be understood without actually solving them. 
Eq. (|33p can be thought as an equation of a point particle 
in one dimension pi l33l| . 



dV 
da 



(36) 
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moving in potential field 

{ 4irG Pb a 2 Aa 2 \ 
V(a) = -\— r - + —), (37) 

where we have used the fact that pb ~ a~ 3 . The Hub- 
ble equation acquires the form of the total energy of the 
mechanical particle 

E = y + V(a) (38) 

where E = —K/2. In order to make the mechanical 
analogy transparent, let us compute the minimum of the 
kinetic energy. If the minimum exists, it should obviously 
correspond to the numerical value of the scale factor that 
gives rise to the maximum of the effective potential V{a). 
It is easy to see that the kinetic energy is minimum if 

a — O'm-! 

a m = {A/k) 1 ' 3 (39) 

(y) = \ ( aV3a1/3 k ) ' (40) 

where A — AnGp^ a 3 ,. Note that V(a) is maximum at 
a = a m . From Eq. (|40[) , we infer that the kinetic energy 
of the system at the top of the potential is, 

/d 2 \ K 3 
( T ^>_ 0i fA>A cSw (41) 

In case A = A c , the system barely makes to the hump of 
the potential (d = 0) corresponding to a m = ao where 
a = as it should be (see Eq. (j3l)]) which is nothing 
but the Einstein's static solution. We are now ready 
to provide the qualitative description of solutions of 
evolution equations. For A < A c , the kinetic energy 
is formally negative for a = a m which means that it 
vanishes before the particle reaches the maximum of 
the potential. In Figf2] we have displayed the plot of 
V(a) versus the scale factor a. We show three possible 
configuration of interest: (A) Corresponds to motion 
starting from the left of the barrier with a = 0. (B) 
Depicts the situation in which the potential barrier is 
approached beginning from the right with a large value 
of the scale factor. (C) Represents the possibility of 
static solution. 

We first analyse the case of K > or E < which 
gives rise to a variety of interesting possibilities. 
1. A < A c : In this case, the kinetic energy is insufficient 
to overcome the potential barrier giving rise to the 
following interesting solutions. 

• (a) Oscillating solution: In this case, motion starts 
from a = with insufficient kinetic energy to reach 
the hump of the potential. In this situation, the scale 
factor increases up to a maximum value where d = 
for a < a rn marking the turning point followed by the 
contraction to a = 0. 





• 




V(a) 









C 


> 




Ati 





FIG. 2: Plot of the effective potential V(a) versus the scale 
factor a. Configurations (A) & (B) correspond to motion of 
system beginning from a — and a = oo respectively. (C) 
corresponds to static solution unstable under small fluctua- 
tions. 



• (b) Bouncing universe: If the potential barrier is 
approached from the right side with a = oo, the scale 
factor first decreases and reaches a minimum value and 
then bounces to expanding phase as the kinetic energy 
is not enough to overcome the barrier. 

• (c) Einstein static universe : This configuration 
corresponds to the maximum of the potential with 
d = and d = 0, possible for a particular value of A, 
obtained earlier. Clearly, static universe corresponding 
to point particle sitting on the hump of the potential, 
is not stable. Small perturbations would derive it to ei- 
ther contracting (a — * 0) or expanding (a — * oo) universe. 

2. A > A c : The kinetic energy is sufficient to overcome 
the barrier for this choice of A. As a result, motion first 
decelerates till the system reaches the top of the potential 
and then slides down the hill with acceleration. Scale fac- 
tor exhibits the point of inflection at a(t) = a m < ao. If A 
slightly exceeds its critical value, an interesting possibil- 
ity dubbed loitering universe can be realized. The scale 
first increases as it should, approaches ao and remains 
nearly frozen for a substantial period before entering the 
phase of acceleration. Such a scenario has important im- 
plications for structure formation. 

For K < or E > 0, the system always has enough 
kinetic energy to surmount the barrier allowing the scale 
factor to increase from a = to large values as time 
increases. This case is similar to the one with K > and 
A > A c . 

For any given value of A, the scale factor exhibits the 
point of inflection at a = a m = (AirGp^ a^/ A) 1 / 3 . This 
also clear from Eq. (|3"2")) & (|33|) in which the first term is 
of attractive character and dominates in the beginning 
leading to deceleration. However, as the scale factor in- 
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creases and reaches a particular value, the repulsive term 
takes over; the scale factor exhibits the point of inflection 
and the expansion becomes accelerating thereafter. 

Observations should tell us when deceleration changed 
into acceleration. This crucially depends upon how 
AnGpfr compares with A or how p\ compares with 

pf) /2. The transition from deceleration to acceleration 
should have taken place around the present epoch. Had 
it happened much earlier it would have obstructed struc- 
ture formation 3J] . We shall come back to this point to 
confirm that cosmic acceleration is indeed a recent phe- 
nomenon. 



IV. BEYOND NEWTONIAN PHYSICS: 
PRESSURE CORRECTIONS 

The formalism of Newtonian cosmology is not appli- 
cable to relativistic fluids. Relativistic fluids essentially 
have non-zero pressure. For instance, radiation is a rel- 
ativistic fluid with pressure Pb = pbC 2 /3. The cosmolog- 
ical constant also belongs to the category of relativistic 
systems. In general theory of relativity, pressure appears 
on the same footing as energy density. Here we present 
heuristic arguments to capture the pressure corrections 
in the evolution equations (see Ref . [22| ) . 

Let us consider a unit comoving volume in the expand- 
ing universe and assume the expansion to be adiabatic. 
The first law of thermodynamics then tells that 



dE + P b dV = 



(42) 



where P b (t) is the pressure of background fluid. The 
first law of thermodynamics applies to any system, be 
it relativistic or non-relativistic, classical or quantum — 
thermodynamics is a great science. 

The energy density of the fluid can always be expressed 
through the mass density, 



E 



4tt 



3 2 

a p b c 



(43) 



Substituting (|4"3"|) into (|4"2")l . we obtain the continuity 
equation in the expanding universe, 



p b + 3H p b 



P, 







(44) 



Thus the continuity equation responds to pressure cor- 
rections: pb —* pb + Pb/c 2 . For a non-relativistic fluid, 
rest energy density dominates over pressure and the sec- 
ond term in the parenthesis can be neglected. For in- 
stance, for dust, Pb ~ 0. At early times, universe was 
hot and was dominated by radiation. Hence the early 
universe should be treated by relativistic theory; Newto- 
nian description becomes valid at late times when matter 
dominates. For the sake of convenience, we shall use the 
unit c = 1. With this choice, relativistic mass density 
and energy density are same. 



We can now present cosmological constant as a per- 
fect fluid with constant energy density. The continuity 
Eq. (j44j) then implies that p\ = —Pa- Next, we claim 
that the correct equation of acceleration in case of back- 
ground fluid with energy density pb and pressure Pb is 
given by 



{p b + 3P 6 ; 



A 
3~ 



To verify, let us multiply Eq. (|4"5|) left right by d 



1 d 
2~dt 



AttG . . A 

■— — a (pba + oPbaj + —aa 



(45) 



(46) 



Using the continuity equation, we can express the term 
containing pressure Pb in Eq. (|46[) through pb, p\ and d 



1 d 
Ydt 



(a 2 ) 



4ttG d 
3 ~di 



p b a 



A 
— < 

6 



(47) 



which can put in form of Friedmann equation in the pres- 
ence of matter with non-zero pressure. 



8ttG . , K A 



(48) 



3 r ' ' a 2 3 

We again observe that pressure corrects the energy den- 
sity. Positive pressure adds to deceleration where as 
the negative pressure contributes towards acceleration, 
see Eq.([45|). It looks completely opposite to our intu- 
ition that highly compressed substance explodes out with 
tremendous impact whereas in our case pressure acts in 
the opposite direction. It is important to understand 
that our day today intuition with pressure is related to 
pressure force or pressure gradient. In a homogeneous 
universe pressure gradients can not exist. Pressure is a 
relativistic effect and can only be understood within the 
frame work of general theory of relativity. Pressure gra- 
dient might appear in Newtonian frame work in the in- 
homogeneous universe but pressure can only be induced 
by relativistic effects. Strictly speaking, it should not 
appear in Newtonian cosmology. This applies to A also 
with negative pressure which we introduced in Newto- 
nian cosmology by hand. Eas.([S )) .(|35] l & |gSJ) coincide 
with the evolution equations of relativistic cosmology. 
Their derivation presented here is heuristic. The rigorous 
treatment can only be given in the framework of general 
theory of relativity where cosmological constant appears 
naturally. 

In order solve the evolution equations, we need a rela- 
tion between the energy density and pressure known as 
equation of state. In case of barotropic fluid the equa- 
tion of state is given by Wb = Pb/ Pb- Dust and radiation 
correspond to iu& = 0,1/3 respectively. Assuming that 
universe is filled with perfect fluid with constant equa- 
tion of state parameter w b , we find from Eqs. (|44]) & (|34|) 
in case of K = 0, 



p b oc a - 3 ( 1+ ™) 

a(t) oc f3(i+"<> , [w > -1) 

a(t) cx eV^* ( w = -1) 



(49) 
(50) 
(51) 
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In case of radiation. Wb = 1/3 and as a result pb = p r oc 
a -4 . In contrast to the case of dust dominated universe, 
the radiation energy density decreases faster with the 
expansion of universe. The positive radiation pressure 
adds to energy density making the gravitational attrac- 
tion stronger. Consequently, the Hubble damping in the 
conservation equation increases allowing the energy den- 
sity decrease faster than dust in expanding universe. This 
can also be understood in a slightly different way, if we 
assume that radiation consists of photons. As universe 
expands, the number density of photons scales as a -3 as 
usual. But since any length scale in the expanding uni- 
verse grows proportional to the scale factor, the energy 
of a photon, hc/X decreases as 1/a leading to p r ~ a -4 
and a(t) oc t 1 / 2 . It is clear that radiation dominated at 
early epochs as pm ~ a -3 for dust. 

Let us make an important remark on the dynamics 
in the early universe which was dominated by radiation 
(for simplicity, we ignore here other relativistic degrees 



of freedom). As p r 



, the first term on the RHS of 



evolution of Hubble equation dominates over the curva- 
ture term K/a 2 ; obviously, cosmological constant plays 
no role in the present case. We therefore conclude that 
all the models effectively behave as K = model at early 
times, 



8ttG ( )Oo 
3 Pr a 4 



a(t) 
ao 



32irGp 



(or 



1/4 



,1/2 



(52) 



We next assume that radiation was in thermal equilib- 
rium characterized by the black body distribution, 



(53) 



Pr = bT\ 



where b is the radiation constant. From Eqs. ([52|) & (j53|) . 
we find how temperature scales with the expansion of 



universe, 



T 



1/4 



a 



(54) 



which on using Eq. (|52j) tells us how early universe cooled 
with time. 



3b 



(55) 



At t = 0, both the radiation density and temperature be- 
come infinitely large; all the physical quantities diverged 
at that time referred to as big bang. The big bang singu- 
larity is not the artifact of homogeneity and isotropy. It is 
a generic feature of any cosmological model based upon 
classical general theory of relativity. Classical physics 
breaks down as big bang is approached. In the frame- 
work of classical general relativity, the big bang is taken 
to be the beginning of our universe. Universe was thus 
born in a violent explosion like event throwing away cos- 
mic matter and giving rise to expansion of universe. Since 



gravity is attractive (provided universe is filled with mat- 
ter of non- negative pressure), its roll is to decelerate the 
expansion. What caused big bang, has no satisfactory 
answer. The big bang is a physical singularity which 
should be treated by quantum gravity. The inflationary 
paradigm can mimic big bang without singularity but 
in that case, we do not know what caused inflation! In 
the cosmic history, there was an epoch when matter took 
over leading to matter dominate era. It turns out that 
it took around 10 5 years for radiation energy density to 
equalize with energy density of matter. The age of uni- 
verse, i.e., the time elapsed since the big bang till the 
present epoch given by (|30[) changes insignificantly, if we 
consider the universe filled with radiation and dust both. 
This is because the time taken from the big bang till radi- 
ation matter equality is negligibly small as compared to 
the actual age of universe which is around 14 Gyr. Thus 
the age given by (f3T))) is a reliable theoretical estimate. 
Unfortunately, the age given by Eq. (f3"0|) falls short than 
the age of some very old objects found in the universe. 
This is one of the old problems of hot big bang model. 
We shall discuss its possible remedy in the dark energy 
dominated universe. 



A. Dark energy 

Eqs. (|4"8"|) and (|4"6")) tell us that the positive cosmologi- 
cal constant A contributes positively to the background 
energy density and negatively to pressure. It can be 
thought as a perfect barotropic fluid with, 



PA 



A 

8^G ; 



Pa = - 



A 

8^G 



(56) 



which corresponds to wa = — 1. In general, we find from 
Eq. (|45|) that expansion has the character of acceleration 
for large negative pressure, 



a 
a 



4irG 



(Pb + 3P h ) 



(57) 



a>0=>Pf,<— : Dark energy. 



where we have included A in the background fluid. Thus 
we need an exotic fluid dubbed dark energy to fuel the ac- 
celerated expansion of universe. The various data sets of 
complimentary support the late time acceleration of uni- 
verse. The simplest candidate of dark energy is provided 
by the cosmological constant with w\ = — 1 Observations 
at present do not rule out the phantom dark energy with 
w < — 1 corresponding to super acceleration. In this case 
the expanding solution takes the form, 



o(t) = % - ty 

H 



u-t 



(n = 2/3(1 + w) <0) (58) 
(59) 



where t s is an integration constant. It is easy to see that 
phantom dominated universe will end itself in a singu- 
larity, in future known as big rip or cosmic doomsday as 
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t — > t s . Clearly, as t — > i a , both the Hub ble p arameter 
and the background energy density diverge [36J. 



B. Age crisis and its possible resolution 

Apart from the cosmic acceleration, dark energy has 
important implications, in particular, in relation to the 
age problem. In any cosmological model with normal 
form of matter, the age of universe falls short compared 
to the age of some known objects in the universe. Since 
the age of universe crucially depends upon the expansion 
history, it can serve as an important check on the model 
building in cosmology. In order to appreciate the prob- 
lem, let us first consider the case of flat dust dominated 
Universe (TIm = 1) hi which case as shown earlier, 
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(60) 



FIG. 3: Plot of age of Universe versus Qm (at present epoch) 
for a flat universe with matter and dark energy with constant 
equation of state parameter w, from Ref. [Ic| 



The observational uncertainty of Hq gives rise to the fol- 
lowing estimate, 

Hq 1 = 9.8/i _1 Gj/r (61) 
0.64 < h < 0.8 -> t = (8 - 10)Gyr (62) 



matter dominated era and we, therefore, have omitted fi r 
in Eq. (|63|) . In case dark energy is cosmological constant 
(wa = —1), we get the analytical expression for the age 
of Universe, 



This model is certainly in trouble as its prediction for 
age of Universe fails to meet the constrain following 
from the study of ages of old stars in globular clusters: 
\2Gyr < to < 15Gyr[35||. One could try to address the 

problem by invoking the open model with fljJJ-' < 1. In 
this case the age of universe is expected to be larger 
than the flat dust dominated Universe — for less amount 
of matter, it would take longer for gravitational attrac- 
tion to slow down the expansion rate to its present value. 
Looking at Eq. ((22|) . it is not difficult to guess that in 

this case, H t -> 1 for -> which is a substan- 
tial improvement. However, this model is not viable for 
several reasons. In particular, the study of large scale 
structure and its dynamics constrain the matter density: 
0.2 < fly < 0.3 and observations on CMB un-isotropy 
reveal that universe is critical to a good accuracy. 

The age problem can be resolved in a flat universe 
dominated by dark energy. Let us rewrite the Friedmann 
equation in a convenient form, 



fOo 

V a 



o (0) (ao 

il DE{ a 



3(1+™) 



(63) 



which allows us to write the expression of to in the closed 
form 



to 



dz 



1 



Ho Jo (i + z ) [n$(i + zf + (i + *)3(i-h») 



1/2 



(64) 

where fl^j is the contribution of dark matter and (1 + 
z ) = a o/°5 z being the rcdshift parameter. The domi- 
nant contribution to the age of universe comes from the 



to 



2H n 



"A 



In 



1 + fi 



1/2 



(0)!/2 



(65) 



n 



M 



For dark energy other than the cosmological constant, 
the integral in Eq. (|64p should be computed numerically. 
In FigOU we have plotted the age of universe versus the 
r various possibilities of dark energy including the 
phantom one. The age constraint can be met by flat dark 



energy models provided that —2 



-0.5 for Q 



(0) 
M 



lying between 0.2 and 0.3, see Ref .]l3]. It is remarkable 
that hot big bang model can be rescued by introducing 
the dark energy component. Interestingly, cosmological 
constant was invoked to address the age problem before 
the invention of cosmic acceleration. The observation of 
cosmic acceleration in 1998 was a blessing in disguise for 
cosmological constant. 



The discovery of cosmic acceleration and its 
confirmation 



The direct evidence of current acceleration of universe 
is related to the observation of luminosity distance by 
high redshift supernovae by two groups independently in 
1998. The luminosity distance for critical universe domi- 
nated by non-relativistic fluid and cosmological constant 
is given by 



d L = 



(1 + *) 
Ho 



dz' 



(i + zf + r4°ki + z') 3(1+M,) 

(66) 
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Eq. (|66|) is the expanding universe generalization of abso- 
lute luminosity L s of a source and its flux J 7 at a distance 
d given by T = L s /(Aird 2 ). It follows from Eq. ([66| that 
Dl — z/Hq for small z and that 

d L = 2 (l + z - (1 + 2 ) 1/2 ) # -\ 0$ = 1 (67) 
d L =z{\ + z)H^\ Q^ E = Q A = 1 (68) 

which means that luminosity distance at high rcdshift 
is larger in universe dominated by cosmological constant 
which also holds true in general for an arbitrary equa- 
tion of state w corresponding to dark energy. Therefore 
supcrnovac would appear fainter in case the universe is 
dominated by dark energy. The luminosity distance can 
be used to estimate the apparent magnitude m of the 
source given its absolute magnitude M 

m-M = 51ogf — — ) + 25 (69) 
\ Mpc J 

Let us consider two supernovae 1997ap at redshift z = 
0.83 with m = 24.3 and 1992p at z = 0.026 with M = 
16.08 respectively. Since the supernovae are assumed to 
be the standard candles, they have the same absolute 
magnitude. Eq. (|6"9")) then gives the following estimate 

H d L ~ 1.16 (70) 
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FIG. 4: Plot of the luminosity distance HqcLl versus the red- 
shift z for a flat cosmological model. The black points come 
from the "Gold" data sets by Riess el al. [38|], whereas the 
red points show the recent data from HST. Three curves show 
the theoretical values of HqcIl for (i) Q$ — 0, = 1, (ii) 
= 0.31, Q A = 0.69 and (iii) £1$ = 1, fi A = 0. From 
Ref. HI. 



Then theoretical estimate for the luminosity distance is 
given by 



dL 
d L 



0.95# , 
1.23if -\ 



(0) 
M 

(0) 
M 



0.3, n A = 0.7 



(71) 
(72) 



where we have used the fact that, d^ — z/Hq for small 
z. The above estimate lands a strong support to the 
hypothesis that late time universe is dominated by dark 
energy (see, Fi^4j . 

The observations related to CMB and large scale struc- 
ture (LSS) provide an independent confirmation of dark 
energy scenario. The acoustic peaks of angular power 
spectrum of CMB temperature anisotropics contains im- 
portant information. The location of the major peak tells 
us that universe is critical to a good accuracy which fixes 
for us the cosmic energy budget. Specializing the Fried- 
mann Eq. (|6"3")) to the present epoch (a = an), we have 



O 



(0) 



O 



(o) 

M 



(0) 
DE 



(73) 



The contribution of radiation to total fractional energy 
density OS; i s negligible at present. The study of large 
scale structure and its evolution indicate that nearly 
30% of the total energy content is contributed by non- 
luminous component of non-barionic nature with dust 
like equation of state popularly known as dark matter. 
The missing component which is about 70% is dark en- 
ergy. The recent data on baryon acoustic oscillation is 



yet another independent probe of dark energy. The com- 
bined analysis of data of complimentary nature demon- 



strate that f2 



(0) 

DE 



0.7 and f2 



(0) 
M 



0.3, see FigEl The 



constraint on the equation of state parameter w and 
f2^ shows that w is restricted to a narrow strip around 
w\ = — 1 (Figl/aJ) . It is clear from the figure that the 
combined analysis allows super-negative values of w cor- 
responding to phantom energy. Let us now confirm that 
the transition from deceleration to cosmic acceleration 
took place in the recent past. Indeed, observations allow 
to estimate the time of transition from deceleration to ac- 
celeration. Let us rewrite Eq. (|38[) through dimcnsionlcss 
density parameters, 



Hi 



si M "o 



(74) 



Using Eq (|74p . we can find out the numerical value of 
(a /ao) corresponding to the minimum of kinetic energy 
(d 2 /2) which precisely gives the transition from deceler- 
ation to acceleration, 



ao 



(0)' 



1/3 



2fi/ 



1/3 



2ft A 



Ztr 



O 



(0) 
M . 



1 ~ 0.67 



(75) 



for the observed values of density parameters (O^ ~ 0.3; 
51a — 0.7) and this confirms that the contribution of 
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FIG. 5: Figure shows the best fit regions in the Qm) 
plane obtained using the CMB, BAO and supernovae data, 
from Ref.fH 




FIG. 6: Constraints on the dark energy equation of w and 
Qm obtained from CMB, BAO and supernovae observations, 
from Ref.f37|l 



A to cosmic dynamics became important at late times 
such that the cosmic acceleration is indeed a recent phe- 
nomenon. 



V. RELATIVISTIC COSMOLOGY 

In the last section we presented heuristic arguments to 
capture the pressure effects in the evolution equations. 



Pressure in cosmology is a relativistic effect which can 
be consistently understood in the frame work of general 
theory of relativity. Einstein equations are complicated 
non-linear equations which do admit analytical solutions 
in presence of symmetries. Homogeneity and isotropy 
of universe is an example of a generic symmetry of space 
time. The assumption of homogeneity and isotropy forces 
the metric to assume the FRW form 



ds 2 = -dt 2 + a 2 (t) 
K = 0, ±1 



dr 2 



1-Kr 2 



r 2 (dd 



2 ' sin 2 



(76) 



where a(t) is scale factor. Coordinates (r,0, 4>) are the 
comoving coordinates. A freely moving particle comes to 
rest in these coordinates. 

Eq. (|76[) is purely a kinemetic statement. The infor- 
mation about dynamics is contained in the scale factor 
a(t). Einstein equations allow to determine the scale fac- 
tor provided the matter contents of universe is specified. 
Constant K in the metric (|76|) describes the geometry of 
the spatial section of space time. K = 0, ±1 corresponds 
to spatially flat, sphere like and hyperbolic geometry re- 
spectively. 

The differential equation for the scale factor follows 
from Einstein equations 



8irGT uv 



(77) 



where G^ is the Einstein tensor, i? /t „ is the Ricci tensor. 
The energy momentum tensor T^ v takes a simple form 
reminiscent of ideal perfect fluid in FRW cosmology 



T» = Diag(-p b ,P b ,P b ,P b ) 



(78) 



Note that pressure in general theory of relativity appears 
on the same footing as energy density. In the FRW back- 
ground, the components of G^ v can easily be computed 



Lr 



1 



K), G\ = — (2ad + a 2 + K) (79) 



Other components of G ' ^ v are identically zero. Einstein 
equations then give rise to the following two independent 
equations 



8ttG 
4ttG 



Pb 



K 
772" 



(80) 
(81) 



We remind that p b designates the total energy density 
of all the fluid components present in the universe. The 
continuity equation p b + 3H(p b + P b ) = can be obtained 
by using Eqs. (|80"|) & (f8"Tj) which also follows naturally 
from the Bianchi identity. As mentioned earlier, we can 
normalize the scale factor to a convenient value at the 
present epoch in case of specially flat geometry. In other 
cases, it should be determined from the relation o,qHq = 
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Qfig - 1|) where defines the total energy content of 
universe at the present epoch. 

Let us note that the Einstein equations ([77)) with the 
energy momentum tensor of standard fluid with positive 
pressure can not lead to accelerated expansion. The re- 
pulsive effect can be captured either by supplementing 
the energy momentum tensor (on right hand side of Ein- 
stein equations) with large negative pressure or by mod- 
ifying the geometry itself. i.,c. the left hand side of Ein- 
stein equations. We can ask for a consistent modification 
of Einstein equations (equation of motion should be of 
second order with the highest derivative occurring lin- 
early so that the Cauchy problem is well posed) in four 
space time dimensions within the classical frame work. 
Under the said conditions, the only admissible modifica- 
tion is provided by the cosmological constant. Thus we 
can add a term Kg^ v on the left hand side of Eq. (|77|) 
which we can formally carry to the right hand side and 
interpret it as the part of energy momentum tensor of a 
perfect fluid[I2](see also Refs.fioj for a different approach 
to cosmological constant), 



G,, 



Rl 



-g^R = 8ttGT^ - K 9llv (82) 



Such a modification is allowed by virtue of Bianchi iden- 
tity. It is remarkable that cosmological constant does not 
need adhoc assumption for its introduction; it is always 
present in Einstein equations. It could be considered as a 
fundamental constant of classical general theory of rela- 
tivity at par with Newton's constant G. It is also interest- 
ing to note that model based upon cosmological constant 
is consistent with all the observational findings in cos- 
mology at present. However, there are deep theoretical 
problems related to cosmological constant. 



A. Theoretical issues associated with A 

There are important theoretical issues related to cos- 
mological constant. Cosmological constant can be as- 
sociated with vacuum fluctuations in the quantum field 
theoretic context [HI, HH, 0] • Though the arguments are 
still at the level of numerology but may have far reaching 
consequences. Unlike the classical theory, the cosmologi- 
cal constant in this scheme is no longer a free parameter 
of the theory. Broadly the line of thinking takes the 
following route. The ground state energy dubbed zero 
point energy or vacuum energy p vac of a free quantum 
field with spin j given by 



= ± (-l) 2j (2j + 1) j g/^W (83) 

\ 2 3 



d 3 k 



Air 2 



(84) 



is ultraviolet divergent. This contribution is related the 
ordering ambiguity of fields in the classical Lagrangian 
and disappears when normal ordering is adopted. Since 



this procedure of throwing out the vacuum energy is ad- 
hoc, one might try to cancel it by introducing the counter 
terms. The later, however requires fine tuning and may 
be regarded as unsatisfactory. The divergence is related 
to the modes of very small wavelength. As we are ig- 
norant of physics around the Planck scale, we might be 
tempted to introduce a cut off around the Planck length 
L p and associate with this a fundamental scale. Thus 
we arrive at an estimate of vacuum energy p vac ~ M 4 

(corresponding mass scale- M vac ~ plat) which is away 
by 120 orders of magnitudes from the observed value of 
this quantity which is of the order of 10~ 48 (GeU) 4 . The 
vacuum energy may not be felt in the laboratory but 
plays important role in GR through its contribution to 
the energy momentum tensor as 

< T^v > = -/WflV, Pvac = A/8irG (85) 

and appears on the right hand side of Einstein equations. 

The problem of zero point energy is naturally resolved 
by invoking supersymmetry which has many other re- 
markable features. In the supersymmetric description, 
every bosonic degree of freedom has its Fermi counter 
part which contributes zero point energy with opposite 
sign compared to the bosonic degree of freedom thereby 
doing away with the vacuum energy. It is in this sense the 
supersymmetric theories do not admit a non-zero cosmo- 
logical constant. However, we know that we do not live in 
supersymmetric vacuum state and hence it should be bro- 
ken. For a viable supersymmetric scenario, for instance 
if it is to be relevant to hierarchy problem, the supper- 
symmetry breaking scale should be around M susy ~ 10 3 
GeV. We arc still remain away from the observed value 
by many orders of magnitudes. We do not know how 
Planck scale or SUSY breaking scales is related to the 
observed vacuum scale! 

At present there is no satisfactory solution to cosmo- 
logical constant problem. One might assume that there 
is some way to cancel the vacuum energy. One can then 
treat A as a free parameter of classical gravity similar to 
Newton constant G. However, the small value of cosmo- 
logical constant leads to several puzzles including the fine 
tuning and coincidence problems. The energy density in 
radiation at the Planck scale is of the order of Plank en- 
ergy density pp ~ 10 72 GeU 4 and the observed value of 
the dark energy density, p\ ~ 0.7 x p^ — 10~ 48 GeU 4 
which implies that p\/pp ~ 10~ 120 . Thus p\ needs to 
be fine tuned at the level of one part in 10~ 120 around 
the Plank epoch, in order to match the current universe. 
Such an extreme fine tuning is absolutely unacceptable 
at theoretical grounds. Secondly, the energy density in 
cosmological constant is of the same order as matter en- 
ergy density at the present epoch. The question what 
causes this coincidence has no satisfactory answer. 

Efforts have recently been made to understand A 
within the frame work of string theory using flux com- 
pactification. String theory predicts a very complicated 
landscape of about 10 500 de-Sitter vacua 14j. Using An- 
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FIG. 7: Cosmologically viable evolution of field energy density 
versus the scale factor. The dotted line shows the evolution 
of background (matter/radiation) energy density. The field 
energy density (with different initial conditions) joins the 
scaling regime and mimics the background. At late times 
it exits the scaling matter regime to become the dominant 
component and to account for the late time acceleration. 
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FIG. 8: Evolution of and pb in absence of scaling regime in 
case of overshoot and undershoot. The field remains trapped 
in the locking regime till its energy density becomes compara- 
ble to that of the background component. It then starts evolv- 
ing slowly and overtakes the background to become dominant 
at late times. 

these system viable to cosmology. 



thropic principal, we are led to believe that we live in one 
of these vacua! 

A novel approach to cosmological constant problem is 
provided in Ref.[l5[. The line of thinking takes follow- 
ing route: In the conventional framework, the equations 
of motion for matter fields are invariant under the shift 
of the matter Lagrangian by a constant while gravity 
breaks this symmetry. Thus, one cannot obtain a sat- 
isfactory solution to the cosmological constant problem 
until the gravity is made to respect the same symmetry. 
An effective action suggested by Padmanbhan in Ref . [l5[ 
is explicitly invariant under the " shift symmetry" . In 
his approach, the observed value of the cosmological con- 
stant should arise from the energy fluctuations of degrees 
of freedom located in the boundary of a spacetime region. 



VI. SCALAR FIELD DYNAMICS RELEVANT 
TO COSMOLOGY 

The fine tuning problem associated with cosmo- 
logical constant led to the investigation of cosmo- 
logical dynamics of a variety of scalar field systems 
such as quintessence, phantoms, tachyons and K- 
essencefilH \vj\ , l4ll . l42l [4^ (see review@ for details). 
Scalar fields can easily mimic dark energy at late times 
and posses rich dynamics in the past. We should note 
that scalar fields models do not address the cosmological 
constant problem, they rather provide an alternative 
way to describe dark energy. The underlying dynamics 
of these systems has been studied in great detail in the 
literature. Scalar fields naturally arise in models of high 
energy physics and string theory. It is worthwhile to 
bring out the broad features their dynamics that make 



A. Quintessence 

A standard scalar field (minimally coupled to gravity) 
capable of accounting for the late time cosmic accelera- 
tion is termed as quintessence. Its action is given by 



S = 



(86) 

The energy momentum tensor corresponding to this ac- 
tion is given by 



T 



d^dycf) - g^ u 



g^d^dpef) + 1/(0) 



(87) 



which gives rise the following expression for energy den- 
sity and pressure in FRW background 



The Euler-Lagrangian equation 



(88) 



for the action 



Sd a (j) 54 
in FRW background acquires the form 



(89) 
(90) 



32^+^ = 



(91) 



which is formally equivalent to the continuity equation 
and can put in the form 



5 °exp^|3(l + W (0))^, 



(92) 
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where w(4>) = P$/ ' p$. Eq. (|88|) tells us that for a steep 
potential <\> 2 >> 1^(0), the equation of state parameter 
approaches the stiff matter limit, w{4>) — > 1 where as 
w(4>) — > —1 in case of a flat potential, <f> 2 « V(4>). 
Hence the energy density scales as ~ a~ n , < n < 6. 
Let us note that while the field rolls along the steep part 
of the potential, its energy density p$ scales faster than 

Pr- 

From Eq. (|5Tj) we find that 

a > -> Pb + 3P b < => ft < V(<j)) (93) 

which means that we need nearly flat potential to account 
for accelerated expansion of universe such that 

In case of field domination regime, the two conditions in 
Eq. ([94|) define the slow roll parameters which allow to 
neglect the </> term in equation of motion for </>. In the 
present context, unlike the case of inflation, the evolution 
of field begins in the matter dominated regime and even 
today, the contribution of matter is not negligible. The 
traditional slow roll parameters can not be connected to 
the conditions on slope and curvature of potential which 
essentially requires that Hubble expansion is determined 
by the field energy density alone. Thus the slow roll 
parameters are not that useful in case of late time accel- 
eration, though, Eq. ([M|) can still be helpful. 

The scalar field model aiming to describe dark energy 
should possess important properties allowing it to allevi- 
ate the fine tuning and coincidence problems without in- 
terfering with the thermal history of universe. The nucle- 
osynthesis puts an stringent constraint on any rclativistic 
degree of freedom over and above that of the standard 
model of particle physics. Thus a scalar field has to sat- 
isfy several important constraints if it is to be relevant 
to cosmology. Let us now spell out some of these fea- 
tures in detail, see Ref.0, da] for details. In case the 
scalar field energy density p^ dominates the background 
(radiation/matter) energy p bl the former should rcdshift 
faster than the later allowing radiation domination to 
commence which in tern requires a steep potential. In 
this case, the field energy density overshoots the back- 
ground and becomes subdominant to it. This leads to 
the locking regime for the scalar field which unlocks the 
moment the p^ is comparable to p b . The further course 
of evolution crucially depends upon the form the scalar 
potential. For the non-interference with thermal history, 
we require that the scalar field remains unimportant dur- 
ing radiation and matter dominated eras and emerges out 
from the hiding at late times to account for late time ac- 
celeration. To address the issues related to fine tuning, 
it is important to investigate the cosmological scenarios 
in which the energy density of the scalar field mimics the 
background energy density. The cosmological solution 
which satisfy this condition are known as scaling solu- 



tions, 

& = const. (95) 
Pb 

The steep exponential potential V(<f>) ~ exp(Xcj)/Mp) 
with A 2 > 3(1 + u>b) in the frame work of standard GR 
gives rise to scaling solutions whereas the shallow expo- 
nential potential with A < v2 leads to a field dominated 
solution (il0 = 1). Nucleosynthesis further constraints A. 
The introduction of a new relativistic degree of freedom 
at a given temperature changes the Hubble rate which 
crucially effects the neutron to proton ratio at temper- 
ature of the order of one MeV when weak interactions 
freeze out. This results into a bound on A, namely[||, 

= 3(1 + w b )/\ 2 < 0.13 => A > 4.5. (96) 

In this case, for generic initial conditions, the field ul- 
timately enters into the scaling regime, the attractor of 
the dynamics, and this allows to alleviate the fine tuning 
problem to a considerable extent. The same holds for the 
case of undershoot, see FigJT] 

Scaling solutions, however, are not accelerating as they 
mimic the background (radiation/matter). One therefore 
needs some late time feature in the potential. There are 
several ways of achieving this: (1) The potential that 
mimics a steep exponential at early epochs and reduces 
to power law type V ~ (j> 2p at late times gives rise to 
accelerated expansion for p < 1 /2 as the average equation 
of state < w((f>) >= (p - l)/(p + 1) < -1/3 in this 
caseJU 5H. (ii) The steep inverse power law type of 
potential which becomes shallow at large values of the 
field can support late time acceleration and can mimic 
the background at early time[46|. 

The solutions which exhibit the aforesaid features are 
referred to as tracker solutions. For a viable cosmic evo- 
lution we need a tracker like solution. However, on the 
basis of observations, we can not rule out the non-tracker 
models at present. 

In the second class of models where trackers are absent, 
there are two possibilities. First, if p^ scales faster than 
p b in the beginning, it then overshoot the background 
and enter the locking regime. In case of the undershoot, 
the field is frozen from the beginning due to large Hubble 
damping. In both the cases, for a viable cosmic evolu- 
tion, models parameters are chosen such that ~ p\ 
during the locking regime. Hence at early times, the 
field gets locked (tu(</>) = — 1) and waits for the matter 
energy density to become comparable to field energy den- 
sity which is made to happen at late times. The field then 
begins to evolve towards larger values of w(<j>) starting 
from w(cf>) = — l(see Fig|S|). In this case one requires to 
tune the initial conditions of the field. The two classes of 
scalar fields are called freezing and thawing models @, EH- 
In case of tracker (freezing) models, one needs to tune the 
slope of the field potential. Nevertheless, these are supe- 
rior to thawing models as they are capable addressing 
both the fine tuning and the coincidence problems. 
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Before we proceed further, we should make an honest 
remark about scalar field models in general. These mod- 
els lack predictive power: for a give cosmic history, it is 
always possible to construct a field potential that would 
give rise to the desired evolution. Their merits should 
therefore be judged by the generic features which arise 
in them. For instance tracker models deserve attention 
for obvious reasons. Scalar fields inspired by a funda- 
mental theory such as rolling tachyons are certainly of 
interest. 



energy may arise in this case as a transient phenomenon, 
(ii) V{4>) — ► slower then 1/0 2 for cj> — > oo ; these models 
give rise to dark energy as late time attractor. The two 
classes are separated by V{4>) ~ l/</> 2 which is scaling 
potential with w(4>) = const. These models suffer from 
the fine tuning problem; dynamics in this case acquires 
dependence on initial conditions. 

C. Phantom field 



B. Tachyon field as source of dark energy 

Next we shall be interested in the cosmological dy- 
namics of tachyon field which is specified by the Dirac- 
Born-lnfeld (DBI) type of action given by (see Ref. Q and 
references therein), 



S = j -VWy/T 



4>dn 4>^J—gd A x 



(97) 



where on phenomcnological grounds, we shall consider 
a wider class of potentials satisfying the restriction that 
V(4>) — > as <fi — > oo. In FRW background, the pressure 
and energy density of <f> are given by 



p = -y(0Wi-0 2 



(98) 



(99) 



The equation of motion which follows from (|97|) is 

+ 3#0(l-0 2 ) + y(l-0 2 ) = (100) 

where H is the Hubble parameter 

1 



H2 = ZMI {P * + Pb) 



(101) 



Tachyon dynamics is very different from that of the 
quintessence. Irrespective of the form of its potential 

w{<f>) = 4> 2 - 1 => - 1 < w{4>) < (102) 

The investigations of cosmological dynamics shows that 
in case of tachyon field, there exists no solution which 
can mimic scaling matter/radiation regime. These mod- 
els necessarily belong to the class of thawing models. 
Tachyon models do admit scaling solution in presence of 
a hypothetical barotropic fluid with negative equation of 
state. Tachyon fields can be classified by the asymptotic 
behavior of their potentials for large values of the field: 
(i) V((f>) — > faster then l/(j> 2 for <\> — > oo. In this case 
dark matter like solution is a late time attractor. Dark 



The scalar field models discussed above lead to w{4>) > 
— 1 and can not give rise to super acceleration correspond- 
ing to phantom dark energy with w(<p) < — 1 permitted 
by observations, sec Fi^jfJJ The simplest possibility of 
getting phantom energy is provided by a scalar field with 
negative kinetic energy. Phantom field is nothing but 
the Hoyle-Narlikar's creation field (C-field) which was 
introduced in the steady state theory to reconcile the 
model with the perfect cosmological principle. Though 
the quantum theory of phantom fields is problematic, it 
is nevertheless interesting to examine the cosmological 
consequences of these fields at classical level. Phantom 
field is described by the following action 



S = 



-gd x 



(103) 



Its corresponding equation of state parameter is given by 



1 J.2 



4> 2 + V(4>) 



4> 2 - V(4>) 



(104) 



which tells us the w{4>) < -1 for <p 2 /2 < V(<j>). An 
unusual equation of motion for <f> follows from (|103p 



• dV 

ZHcj> - — = 



(105) 



It should be noted that evolution equation of phantom 
field is same as that of the ordinary scalar field but with 
inverted potential allowing the field with zero kinetic en- 
ergy to rise up the hill. As mentioned earlier, phantom 
energy is plagued with 6117 rip singularity which is charac- 
terized by divergence of the Hubble parameter and cur- 
vature of space time after a finite interval of time. In 
such a situation, quantum effects become important and 
one should include higher curvature corrections to gen- 
eral theory of relativity which can crucially modify the 
structure of the singularity. To the best of our knowl- 
edge, the big rip singularity can be fullyresolved in the 
frame work of loop quantum cosmoloau\4l\ . Big rip can 
also be avoided at the classical level in a particular class 
of models in which potential has maximum. In this case, 
the field rises to the maximum of the potential and ulti- 
mately settle on the top of the potential to give rise to 
de-Sitter like behavior. 

For a viable cosmic history, the phantom energy den- 
sity similar to the case of rolling tachyon should be sub- 
dominant at early epochs. The field then remains frozen 
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till late times before its energy density becomes compara- 
ble to matter energy density. Its evolution begins there- 
after. Clearly, dark energy models based upon phantom 
fields belong to the category of thawing models. 



D. Late time evolution of dark energy 

In the preceding subsections, we have described the 
cosmological dynamics of quintessence, phantoms and 
rolling tachyon. These scalar field models fall into two 
broad categories: (i) Tracker or freezing models in which 
the field rolls fast at early stages such that it mimics 
the background with Wf, = 0. At late times, w{<j>) starts 
deviating from dust like behavior and becomes negative 
moving towards de-Sitter phase as the field rolls down its 
potential, (ii) Non-tracker or thawing models are those 
in which the field is trapped in the locking regime due to 
large Hubble damping such that w(4>) = — 1. And only 
at late times, as becomes comparable to the back- 
ground energy density, the field begins to evolve towards 
larger values of w (</>). As demonstrated by Caldwell and 
Linder (42j, these models occupy narrow regions in the 
(w' = dw/dhi(a),w) plane, 

3w(l + w) < w < 0.2w(l + w) Freezing models. 
1 + w < w < 3(1 + w) Thawing models. 

where the upper and the lower bounds are obtained using 
analytical arguments and numerical analysis of generic 
models belonging to both the classes of models. As 
pointed out earlier (see Fig|5]), combined analysis of dif- 
ferent observations reveal that dark energy equation of 
state parameter lies in the narrow strip around w& = — 1. 
The observational resolution between the two classes of 
the model which is of the order of 1 + w is therefore a 
challenge to future observations. 

As mentioned earlier, the phantom and the tachyon 
dark energy models belong to the class of thawing mod- 
els. In this case, we can simplify the dynamics around the 
present epoch by using the approximation that |l+iy| << 
1 and that the slope of the potential is small. The validity 
of the second approximation can be verified numerically 
in each case. In this scheme of a plausible approxima- 
tion, one arrives at an amazing result: All the differ- 
ent dynamical systems, thawing quintessence, phantom, 
tachyon and phantom tachyon follow a unique evolution- 
ary track. The distinction between the four classes of 
scalar field systems and the distinction between differ- 
ent models within each class is an effect of higher order 
than |1 + w|[43| which certainly throws a great challenge 
to future generation experiments! Indeed, a recent ex- 
amination of observational data including 397 Type la 
supernovae at redshifts 0.015 < z < 1.55 has shown that 
evolving dark energy models provide a slightly better fit 
to the data than the cosmological constant [48| . If future 
data confirms this result then it could mean that cosmic 
acceleration is currently slowing down which may have 
important consequences for dark energy model building. 



E. Quintessential inflation on brane: A beautiful 
model that does not work 



Quintessential inflation refers to attempts to describe 
inflation and dark energy with a single scalar field. The 
unifications of the two phases of accelerated expansion 
could be realized in the framework of Randall-Sundrum 
(RS) brane worlds [45l. [46|. In order to achieve this, the 
field potential should be flat during inflation but steep in 
radiation and matter dominated eras such that p§ could 
mimic the background energy density at early epochs. At 
late times, it should become flat so as to allow the cur- 
rent acceleration of universe. Since the potential does not 
exhibit minimum, the conventional reheating mechanism 
does not work in this scenario. One could employ al- 
ternative mechanisms such as reheating via gravitational 
particle production or instant preheating. It is not real- 
istic to have a potential which changes from flat to steep 
and back to flat at late times(sec, FiglS]). However, it is 
generic to have a potential which is steep and allows to 
track the background at early epochs and gives rise to a 
viable late time cosmic evolution. 

In case of a steep potential, the field energy density 
scales faster that radiation energy density leading to the 
commencement of radiative regime. But a steep poten- 
tial can not support inflation in FRW cosmology. This 
is precisely where the brane assisted inflation comes to 
our rescue. In RS brane world model, the Fricdmann 
equation is modified to, 



8irG I Pb 



(106) 



where As is the brane tension. The presence of quadratic 
density term in the Friedmann equation changes the dy- 
namics at early epochs in crucial manner. Consequently, 
the field experiences greater damping and rolls down its 
potential slower than it would during the conventional in- 
flation. This effect is reflected in the slow-roll parameters 
which have the form, 



£frw - 



1 + V/Xi 
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(107) 
(108) 



where efrw and t/frw are the standard slow-roll pa- 
rameters in absence of brane corrections. The influence 
of brane corrections becomes specially important when 
V/Xb >> 1. In this case, we have, 

€~e FRW {V/X B )'\ t 1 ~2i 1 frw{VX b )- 1 (109) 

which tells us that slow-roll (e, r/ << l)is possible when 
V/Xb » 1 even if the potential is steep {cfrw, Vfrw > 
1). As the field rolls down its potential, the high energy 
brane correction to Fricdmann equation disappears giv- 
ing rise to the natural exit from inflation. 

It is possible to choose potentials suitable to 
quintessential inflation and fine tune the model pa- 
rameters such that the model respects nucleosynthe- 
sis constraints and leads to observed late time cosmic 
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acceleration[45|,|46|. However, the problem occurs on the 
other side. Recent measurements of CMB anisotropics 
place fairly strong constraints on inflationary models. 
The tensor to scalar ratio of perturbations turns out to 
be lager than its observed value in case of steep brane 
world inflation. Clearly, the brane world unification of 
inflation and dark energy is ruled out by observation. 



V(0) 




FIG. 9: A desired form of potential for quintessential infla- 
tion. It is generic to have a steep potential at early times with 
brane corrections helping the slow-roll of the field. 




R-(l/2)g^d^d„<j) 



V(d>) - f(cj))R 2 



where R^ B is the Gauss-Bonnet term, 



RGB = R 2 — ^R l ivR ilV + RafinvR 



(110) 
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The dilaton potential V(<p) and its coupling to curvature 
/(</>) are given by, 



V{4>) 



f(4>) 



>0) 
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The cosmological dynamics of system JllOp in FRW back- 
ground was investigated in Ref. [5l|, |52|. It was demon- 
strated that scaling solution can be obtained in this case 
provided that [i = a. In case (j ^ a, the de-Sitter so- 
lution is a late time attractor. Hence, the string cur- 
vature corrections under consideration can give rise to 
late time transition from matter scaling regime. Un- 
fortunately, it is difficult to reconcile this model with 
nucleosynthesis p3ll |52| constraint. 



DGP model 



VII. MODIFIED THEORIES OF GRAVITY AND 
LATE TIME ACCELERATION 

The second approach to late time acceleration is re- 
lated to the modification of left hand side of Einstein 
equations or the geometry of space time. It is perfectly 
legitimate to investigate the possibility of late time accel- 
eration due to modification of Einstein-Hilbert action In 
the past few years, several schemes of large scale modifi- 
cations have been actively investigated. Some of these 
modifications are inspired by fundamental theories of 
high energy physics where as the others are based upon 
phenomcnological considerations. In what follows, we 
shall briefly describe the modified theories of gravity and 
their relevance to cosmology. 



A. String curvature corrections 

It is interesting to investigate the string curvature cor- 
rections to Einstein gravity amongst which the Gauss- 
Bonnet correction enjoys special status[5(| HH, HI, HI, 
l54l l55l I5H IoTI . HI]. These models, however, suffer from 
several problems. Most of these models do not include 
tracker like solution and those which do are heavily con- 
strained by the thermal history of universe. For instance, 
the Gauss-Bonnet gravity with dynamical dilaton might 
cause transition from matter scaling regime to late time 
acceleration allowing to alleviate the fine tuning and co- 
incidence problems. Let us consider the low energy effec- 



In DGP model, gravity behaves as four dimensional 
at small distances but manifests its higher dimensional 
effects at large distances. The modified Friedmann equa- 
tions on the brane lead to late time acceleration. The 
model has serious theoretical problems related to ghost 
modes and super luminal fluctuations. The combined 
observations on background dynamics and large angle 
anisotropics reveal that the model performs much worse 
than ACDM[6(j. However, generalized versions of DGP 
can be ghost free and can give rise to transient accelera- 
tion as well as a phantom phase [49|. 



C. f(R) theories of gravity 

On purely phenomenological grounds, one could seek 
a modification of Einstein gravity by replacing the Ricci 
scalar in Einstein-Hilbert action by f(R). The action of 
f(R) gravity is given by [13, 



S = 
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The modified Einstein equations which follow from ()113[) 
have the form, 

f'R^u - V„V„/' + (nf - 9fiv = SttGT^. (114) 

which arc of fourth order for a non- linear function f(R). 
Here prime denotes the derivatives with respect to R. 
The Ricci scalar in FRW background is given by 
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(115) 
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which tells us that the modified Eq. (|114p contains de- 
Sitter space time as a vacuum solution provided that 
/(4A) = 2A/'(4A). The f(R) theories of gravity may in- 
deed provide an alternative to dark energy To see this, 
let us write the evolution equations which follow from 
(|114|) in a convenient form 



H z = 
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where pr and Pr are energy density and pressure con- 
tributed by curvature modification 

PR = ^ r ^-'3HRf" (118) 

P R = 2HRf" + Rf" + -(/ - f'R) + f'R 2 (119) 

pa and Pr identically vanish in case of Einstein-Hilbcrt 
action, f(R) = R as it should be. As an example of f(R) 
model let us consider, f(R) = R — a n /R n , where a n is 
constant for given n. In case of a power law solution 
a(t) ~ t n , the effective equation of state parameter can 
be computed as 



1 2 H 
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3(2n + l)(n + l) 



(120) 



Choosing a particular value of n, we can produce a de- 
sired equation of state parameter for dark energy. 

The functional form of f(R) should satisfy certain re- 
quirements for the consistency of the modified theory 
of gravity. The stability of f(R) theory would be en- 
sured provided that, f'(R) >0 and f (R) > which 
means that graviton is not ghost and scalar degree is 
not tachyon. We can understand the stability conditions 
heuristically without entering into their detailed investi- 
gations. From evolution equations (| 1 16[) & (|116p . we see 
that the effective gravitational constant G e // = G/f 
which should be positive or /' > in order to avoid the 
pathological situation. As for the second condition, V. 
Faraoni has given an interesting interpretation 61| : let 
us consider the opposite case when /" < which means 
that G' eff = -f"G/f' 2 > 0. This implies that gravita- 
tional constant increases for increasing value of R mak- 
ing the gravity stronger. In view of Einstein equations, it 
leads to yet larger value of curvature and so on which ul- 
timately leads to a catastrophic situation. Thus we need 
/" to be positive to avoid the catastrophe. 

Let us note that f(R) gravity theories apart from a 
spin two object necessarily contain a scalar degree of 
freedom. Taking trace of Eq. (|114p gives the evolution 
equation for the scalar degree of freedom, 



I (2/' - f'R) 



8irG, 



-T. 



(121) 



It should be noticed that Eq. (|121[) reduces to an alge- 
braic relation in case of Einstein gravity; in general /' 



has dynamics. It is convenient to define scalar function 

4> as, 



= f-l, 



(122) 



which is expressed through Ricci scalar once f(R) is spec- 
ified. We can write the trace equation fEq. (|121[ )) in the 
terms of V and T as 



U(j) 



dV 



IttG 



T. 



(123) 



which is a Klein- Gordon equation in presence of a deriv- 
ing term. Thus (f> is indeed a scalar degree of freedom 
which controls the curvature of space time. 

The effective potential can be evaluated using the fol- 
lowing relation 



dV_ 
dR 



dV 



1 



dR 



(2/ -/'#)/". 



(124) 



Models which satisfy the stability conditions belong to 
two categories: (1) Either they are not distinguishable 
from AC DM or are not viable cosmologically. (ii) Models 
with disappearing cosmological constant: In these mod- 
els, f{R) — ► for R — > and they give rise to cosmologi- 
cal constant in regions of high density and differ from the 
latter otherwise. In principal, these models can be dis- 
tinguished from cosmological constant. Models belong- 
ing to the second category were proposed by Hu-Sawicki 
and Starobinsky [5^, [6^ (see also Ref.[63| on the simi- 
lar theme). The functional form of f(R) in Starobinsky 
paramctrization is given by, 



f{R) = R + XR Q 



1 R 2 

1 + M 



(125) 



Here n and A are positive. And Rq is of the order of 
presently observed cosmological constant, A = 8irGp va ,c- 
The model satisfies the stability conditions quoted above. 

In the Starobinsky model, the scalar field </>, in the 
absence of matter, is given by 



<KR) = - 



2n\R 



71+1 



(126) 



Notice that R — > oo for <j> — ► 0. For a viable late time cos- 
mology, the field should be evolving near the minimum 
of the effective potential. The finite time singularity in- 
herent in the class of models under consideration severely 
constrains dynamics of the field. 



The curvature singularity and fine tuning of 
parameters 

The effective potential has minimum which depends 
upon n and A. For generic values of the parameters, the 
minimum of the potential is close to <j> = 0(see Fig fTU)) 
corresponding to infinitely large curvature. Thus while 
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FIG. 10: Plot of effective potential for n = 2 and A = 1.2. 
The red spot marks the initial condition for evolution. 



the field is evolving towards minimum, it can easily oscil- 
late to a singular point [64Llo5j]. However, depending upon 
the values of parameters, we can choose a finite range of 
initial conditions for which scalar field <fi can evolve to the 
minimum of the potential without hitting the singular- 
ity. We find that the range of initial conditions allowed 
for the evolution of 4> to the minimum without hitting 
singularity shrinks as the numerical values of parameters 
n and A increase. In the presence of matter, the mini- 
mum of the effective potential moves towards the origin. 
In case of the compact objects such as neutron stars, the 
minimum is extremely near the origin and the singularity 
problem becomes really acute (6oll66|. 



Avoiding singularity with higher curvature 
corrections 

We know that in case of large curvature, the quan- 
tum effects become important leading to higher curva- 
ture corrections. Keeping this in mind, let us consider 
the modification of Starobinsky's model [67l. loll. l69j . 



f(R)=R+—R 2 + R Q \ 

Kq 



then 6 becomes 



-1 



1 



, (127) 



2a 



2n\ 



(1 + 1) 



n+l 



(128) 



In case \R\ is large, the first term which comes from 
aR 2 dominates. In this case, the curvature singular- 
ity, R = ±oo corresponds to (f> = ±oo. Hence, in this 
modification, the minimum of the effective potential is 
separated from the curvature singularity by the infinite 
distance in the <p, V((f>) plane. Though the introduction 
of R 2 term formally allows to avoid the singularity but 
can not alleviate the fine tuning problem as the min- 
imum of the effective potential should be near the in 
generic cases. As for the compact objects, Langlois and 
Babichev[7(|(see, Ref.[7l| also on the similar theme) have 
argued that neutron stars can be rescued from singu- 
larity if a realistic equation of state for these objects is 
used though the numerical simulation is yet challenging 



for densities of the order of nuclear matter density. The 
problem deserves further investigation. 

In scenarios of large scale modification of gravity, one 
should worry about the local gravity constraints. The 
f(R) theories are related to the class of scalar tensor 
theories corresponding the Brans-Dicke parameter ui = 
or the PPN parameter 7 = (1 + uj)/(2 + w) = 1/2 
unlike GR where 7 = 1 consistent with observation 
(I7 — 1| < 2.3 x 10~ 5 ). This conclusion can be escaped 
by invoking the so called chameleon mechanism . In 
case, the scalar degree of freedom is coupled to matter, 
the effective mass of the field depends upon the matter 
density which can allow to avoid the conflict with solar 
physics constraints. However, the problem of singularity 
in these models is genuine and should be addressed. 



VIII. SUMMARY 

We have given a pedagogical exposition of physics of 
late time cosmic acceleration. Most of the part of the 
review should be accessible to a graduate student. The 
discussion of Newtonian cosmology is comprehensive and 
reviews the efforts to put the formalism of Newtonian 
cosmology on rigorous foundations in its domain of va- 
lidity. Heuristic discussion on the introduction of cos- 
mological constant and pressure corrections in evolution 
equations is included. The underlying idea leading to 
late time cosmic acceleration is explained without the 
use of general theory of relativity. The basic features of 
cosmological dynamics in presence of cosmological con- 
stant is presented in a simple and elegant fashion mak- 
ing it accessible to non-experts. The review also gives 
the glimpses of relativistic cosmology, contains important 
notes on the dynamics of dark energy and discusses un- 
derlying features of cosmological dynamics of a variety of 
scalar fields including quintessence, rolling tachyon and 
phantom. Special emphasis is put on the cosmic viabil- 
ity of these models; the cosmological relevance of scaling 
solutions is briefly explained. The review ends with a 
discussion on modified theories of gravity as possible al- 
ternatives to dark energy. The treatment is simple but 
conveys the successes and problems of cosmology in the 
frame work of modified theories of gravity. Basic fea- 
tures of f(R) cosmology arc explained avoiding the cum- 
bersome mathematical expressions. The latest develop- 
ments of f(R) theories with disappearing cosmological 
constant are highlighted. The problem of singularities in 
these models and their possible resolution are discussed. 
I hope the review would be helpful to beginners and will 
also be of interest to experts. 
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